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Abstract In the present article we derive an explicit expression for the trun- 
cated mean and variance for the multivariate normal distribution with ar- 
bitrary rect angu l ar do uble truncation. We use the moment generating ap- 
proach of iTallisI l|l96l|) and extend it to general fi, L and all combinations 

Oof truncation. As part of the solution we also give a formula for the bivari- 
ate marginal density of truncated multinormal variates. We also prove an 
| invariance property of some elements of the inverse covariance after trunca- 

tion. Computer algorithms for computing the truncated mean, variance and 
the bivariate marginal probabilities for doubly truncated multivariate normal 
variates have been written in R and are presented along with three examples. 
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1 Introduction 



The multivariate normal distribution arises frequently and has a wide range 
of applications in fields like multivariate regression, Bayesian statistics or the 
analysis of Brownian motion. One motivation to deal with moments of the 
truncated multivariate normal distribution comes from the analysis of special 
financial derivatives ("auto-callables" or "Expresszertifikate") in Germany 
These products can expire early depending on some restrictions of the un- 
derlying trajectory if the underlying is above or below certain call levels. In 
the framework of Brownian motion the finite-dimensional distributions for 
log returns at any d points in time are multivariate normal. When some of the 
multinormal variates X = {x\,. . - ,x/)' ~ N(fl,T,) are subject to inequality 
constraints (e.g. fl, < x, < £>,), this results in truncated multivariate normal 
distributions. 

Several types of truncations and their moment calculation have be en de- 
scrib ed so far, for example the one-sided rectangular truncation x > a (Tallis 



1961), the rather unusual elliptical and ra dial truncati ons a < x'Rx < b (Tallis 



1963) and the plane truncation Cx > p ( Tallislll965|) . Linear constraints like 



a < Cx < b can often be reduced to rectangular truncation by transformation 
of the variables (in case of a full rank matrix C : a* = C a < x < C b = 
b*), which makes the double rectangular truncation a < x < b especially 
important. 

The existing works on moment calculations differ in the number of vari- 
ables they consider (univariate, bivariate, multivariate) and the types of rect- 
angular truncation they allow (single vs. double truncation). Single or one- 
sided truncation can be either from above (x < a) or below (x > a), but 
only on one side for all variables, whereas double truncation a < x < b can 
have both lower and upper truncations points. Other distinguishing features 
of previous works are further limitations or restrictions they impose on the 
type of distribution (e.g. zero mean) and the methods they use to derive the 
results (e.g. direct integration or moment-generating function). Next, we will 
brief ly outline the line of research. 

iRosenbauml jl96l|) gave an explicit formula for the moments of the bivari- 
ate case with single truncation from below in both variables by direct inte- 
gr ation. His results f or the bivariate normal dist r ibutio n ha ve been extend ed 
by lShah and Parikhl (|l964|) . lRegier and Hamdanl (|l97l|) and lMuthen] Jl990|) to 
double truncation. 

For the multivariate case, iTallid d!96lh derived an explicit expression for 
the first two moments in case of a singly truncated multivariate normal den- 
sity with zero mea n vector and the corr elation m atrix R using th e moment 
generating function. lAmemiyal l|l974|) andlLed l|l979|) extended the lTallisI l|l96l|) 
derivation to a general covariance m atrix L and also eyaluat ed t he re l ation- 
ship between the first two moments. iGupta and Tracyl l ll976|) and lLed jl983|) 
gave very simple recursive relationships between moments of any order for 
the doubly truncated case. But since except for the mean there are fewer equa- 
tions than parameters, these recurrent conditions do not uniquely identify 
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Table 1 Survey of previous works on the moments for the truncated multivariate normal distri- 
bution 



Author 



#Variates 



Truncation Focus 



iRosenbauml 11961ft 



lTallislll96lft 

Shah and Parikhl 11964ft 
Regier and Ham danl 



11971ft 
lAmemivalll974ft 



Gupta and Tracvl 11976ft 



Lee 



Lee 



1979) 



1983) 



Leopard and Tallis 
11989ft 

lMuthenl<1990h 



Manjunath/Wilhelm 



bivariate 



bivariate 
bivariate 



multivariate 
multivariate 
multivariate 
multivariate 

bivariate 



single 



multivariate single 



double 
double 



multivariate single 



double 
single 
double 
single 

double 



multivariate double 



moments for bivariate normal variates 
with single truncation, b\ < y\ < 

oo, &2 < J/2 < 00 

moments for multivariate normal vari- 
ates with single truncation from below 
recurrence relations between moments 
an explicit formula only for the case of 
truncation from below at the same point 
in both variables 

relationship between first and second 
moments 

recurrence relations between moments 
recurrence relations between moments 
recurrence relations between moments 
moments for multivariate normal distri- 
bution with single truncation 
moments for bivariate normal distribu- 
tion with double truncation, b\ < j/j < 
fl l/^2 < 1/2 < a 2 

moments for multivariate normal distri- 
bution with double truncation in all vari- 
ables a < x < b 



moments of order > 2 and are therefore not sufficient for the computation of 
the variance and other higher order moments. 

Table Q] summarizes our survey of existing publications dealing with the 
computation of truncated moments and their limitations. Even though the 
rectangular truncation a < x < b can be found in many situations, no ex- 
plicit moment formulas for the truncated mean and variance in the general 
multivariate case of double truncation from below and /or above have been 
presented so far in the literature and are readily apparent. The contribution 
of this paper is to derive these formulas for the first two truncated moments 
and to ex tend and generalize exi s ting results on moment calculations from 
especiallySnii ( |l96l|) :lLe3 l)l983|) ; lLeppard and Tallis! l|l989fr ; lMuthinl |)i990l) . 



The remainder of this paper is organized as follows. Section 2 presents 
the moment generating function (m.g.f) for the doubly truncated multivariate 
normal case. In Section 3 we derive the first and second moments by differen- 
tiating the m.g.f. These results are completed in Section 4 by giving a formula 
for computing the bivariate marginal density. In Section 5 we present two nu- 
merical examples and compare our results with simulation results. Section 6 
links our results to the theory of graphical models and derives some proper- 
ties of the inverse covariance matrix. Finally Section 7 summarizes our results 
and gives an outlook for further research. 
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2 Moment Generating Function 

The d-dimensional normal density with location parameter vector fi G !R d 
and non-singular covariance matrix L is given by 



( 27r )rf/2|2:|i/2 



exp 



(x-fi)'!:- 1 {x-ji)\, x G M fl . (1) 



The pertaining distribution function is denoted by <$> t , r z(x). Correspondingly, 
the multivariate truncated normal density truncated at a and b, in ]R rf , is de- 
fined as 



<?V,e(x) 



-, for a < x < b, 



p{a<X<b}' 

0, otherwise. 



(2) 



Denote a = P{a<X<b}as the fraction after truncation. 

The moment generating function (m.g.f) of a d-dimensional truncated 
random variable X, truncated at a and b, in R d , having the density f(x) is 
defined as the d-fold integral of the form 

m(t) = E (e t>X ) = [ e t>x f(x)dx. 



Therefore, the m.g.f for the density in © is 



m(t) 



cc{2ny 2 \L\V z Ja 



exp 



{x - }i)' ^ l {x - }i) - 2t' x }dx.(3) 



In the following, the moments are first derived for the special case }i = 0. 
Later, the results will be generalized to all ji by applying a location transfor- 
mation. 

Now, consider only the exponent term in 10 for the case }i — 0. Then we 
have 



which can also be written as 



2 2 



x'E -1 * - 2t'x 



(x-$yL-'(x-Z) 



where £ = £f. 



Consequently, the m.g.f of the rectangularly doubly truncated multivari- 
ate normal is 



m(t) = 



exp 



(*-£)' \dx, (4) 
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where T = $t"Lt. 

The above equation can be further reduced to 

For notational convenience, we write equation as 

m{t) = e T <t> aZ (6) 

where 

1 /•&-? f 1 , 1 

E a(27T) rf / 2 |L| 1 / 2 Ja-l \ 2 J 



3 First And Second Moment Calculation 

In this section we derive the first and second moments of the rectangularly 
doubly truncated multivariate normal density by differentiating the m.g.f.. 

Consequently by taking the partial derivative of I0 with respect to f , we 
have 

dm(t) r 9$Vs de T 
In the above equation the only essential terms which will be simplified 

are 

de T T * 
OTl k=i 

and 



dU dt 



1 fi-i Ja. Ja 



f KX (x)dx d ...dxi, (8) 



if 



where a* = a { - ELi <r#h and h * = b i ~ E*=i "U**- 
Subsequently |[§) is 



where 

F,(x) 

J ... ' 1 1+1 ... / d (p aT ,(x 1 ,..,x i _ 1/ x,x i+1/ ..x d )dx d ...dx i+1 dx i _ 1 ...dx 1 . 

J at J a*- , Jar, , J a*. 



d 

(10) 
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Note that at t k = 0, for all k = 1,2, ...,d, we have a* = a { and b* = bj. 
Therefore, Fj(x) will be the z-th marginal density. An especially c onvenient 
way o f computing these one-dimensional marginals is given in Cartinhour 
dl99fl . 

From 10 - © for k = 1,2, d all t k = 0. Hence, the first moment is 

E(Xf) = -^\ t =Q = E (**(«*) - h(h)) ■ (11) 

Now, by taking the partial derivative of with respect to 1 1, we have 

d 2 m(t) _ T 3 2 d> ffi£ | 93> gE 3e r | | ^ T a<D ff£ 

dtjdtj dtjdtj dtj dtj dtjdti dtj dtj 

The essential terms for simplification are 

9 2 e r _ 

and clearly, the partial derivative of (0 with respect to t; gives 



In the above equation merely consider the partial derivative of the marginal 
density Fj-(fl^) with respect to f;. With further simplification it reduces to 







rK+i 






M-i ■ 


> a t+i 


Ja* 



<J k,q°'j,k 

cm 



where 



'i+l / / 9+1 / d 



fc+1 if— 1 £/+l 



and the short form denotes the vector [x\, ..,x k _i,x k+ i, ..x^) 1 in (d — 1)- 
dimensions and x_ k _„ denotes the (d — 2) -dimensional vector (x\, x k _\, 

Xg -i, Xg + i, x^)' for k ^ q. The above equation (Ml is deduced from 
iLed {l979i pp. 167. Note that for all t k - the term F kq {x,y) will be the bi- 
variate marginal density for which we will give a formula in the next section. 
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Subsequently, dFk d [ b ^ can be obtained by substituting a| by h\ . From O - 
dTHI l at all = 0, k = 1, 2, d, the second moment is 

=t ^ + £^ ^ 

" ' r )(,,/< r i,A 



' £ ^ I] ( ^ ) (fit,,? ) - Ffc /I? («*:/ ) 

fc=l (j^fc V a Kk / L v 



(' 



h,q{h,a q )-F k ,q{h,\)) ■ (16) 



Having derived expressions for the first and second moments for double 
truncation in case of fi = 0, we will now generalize to all ji: 
if Y ~ N(ft,Z) with «*<i/< b*, then X = Y - ji ~ N(0,E) with a = 
a* -/i < x < b* = and E(Y) = E(X) + /< and Coz;(Y) = Coi>(X). 
Equations i fTTI l and | [T6ll can then be used to compute £(X) and Cov(X). Hence, 
for general ji, the first moment is 

d 

E(Yi) = £ <r i>k (F k (a k ) - F k (b k )) + n . (17) 
k=l 

The covariance matrix 

Cov(Yj,Yj) = Cov(Xi,Xj) = E(XiXj) - E(X ! -)E(X / ) (18) 

is invariant to the shift in location. 

The equations (flTb and (f]~8l > in combination with <fTTb and (Tl6ll form our 
desired result and allow the calculation of the truncated mean and truncated 
variance for general double truncation. A formula for the term F k ^(x k , Xq), 
the bivariate marginal density, will be given in the next section. 
We have implemented the moment calculation for mean vector mean, covari- 
ance matrix sigma and truncation vectors lower and upper as a function 

mtmvnorm(mean, sigma, lower, upper) 

in the R package tmvtnorm jWilhelm and Maniun athl2010al: ?), where the code 
is open source. In Section 5 we show a usage example for this function. 



4 Bivariate Marginal Density Computation 



In ord er to compu te the bivaria te marginal densit y in this section we mainly 
follow Tallisl < 1961), p. 223 and lLeppard and Tallisl jl989|) who implicitly used 
the bivariate marginal density as part of the moments calculation for single 
truncation, evaluated at the integration bounds. However, we extend it to the 
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doubly truncated case and state the function for all points within the support 
region. 

Without loss of generality we use a z-transformation for all variates x = 
(x\, . . . , x^Y as well as for all lower and upper truncation points a — (a.\, . . . ,a^)' 
and b = (b\, . . . , b d )' , resulting in a N(0, R) distribution with correlation ma- 
trix R for the standardized untruncated variates. In this section we treat all 
variables as if they are z-transformed, leaving the notation unchanged. 

For computing the bivariate marginal density F cj/r (x cj , x r ) with ciq < Xq < 
bq, a r < x r < b r ,q 7^ r, we use the fact that for truncated normal densities the 
conditional densities are truncated normal again. The following relationship 
holds for x s , z s E R d_2 if we condition on Xq = c q and x r = c r (s 7^ q 7^ r): 

«T l (p A (x s ,x q = c q ,x r = c r ;R) = aT 1 q>{c q ,c r ; p qr )q> d ^ 2 {z 5 ;Rqr) , (19) 

where 



Z S = (X S - B sq . r Cq - B sr . q C r ) I y '(I ~ pl q )(l ~ pl r , q ) (20) 

and R qr is the matrix of second-order partial correlation coefficients for s 7^ 
cj 7^ r. f!> sq . r and & S r.q are the partial regression coefficients of x s on x q and x r 
respectively and p sr .q is the partial correlation coefficient between x s and x r 
for fixed x q . 

Integrating out (d — 2) variables x s leads to F q ,r (x q , x r ) as a product of a bivari- 
ate normal density <p(x q , x r ) and a (d — 2)-dimension normal integral 

F (x — c x — c ) — I**- 1 [ bq+1 



f a R{x s ,Cq,C r )dx s 

flr+1 Ja d 

= or l (p{c q ,C r }p v )®cl-2{AU'B q rs,'R q r) (21) 

where A% and B% denote the lower and upper integration bounds of <&d-2 
given x q = c q and x r = c r : 

Al s = fa - B sq ,Cq ~ B sr . q C r )/ ^(1-P%)(l-Pl r .q) (22) 

Bl = (b s - B sq ,c q - B sr . q c r ) I ^(l-p2 q )^- p 2 r q ). (23) 

The computation of Fq ir (xq,x r ) just needs the evaluation of the normal inte- 
gral in d — 2 dimensions, which is readily available in most statistics 
software packages, for example as the function pmvnorm ( ) in the R package 
mvtnorm (?). 

The bivariate marginal density function 

dtmvnorm(x, mean, sigma, lower, upper, margin=c(q,r) ) 

is also part of the R package tmvtnorm (|Wilhelm and Manjunath l2010al ?), 
where readers can find the source code as well as help files and additional 
examples. 
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bivariate marginal density (x^) 




-2 



Fig. 1 Contour plot for the bivariate truncated density function 



Marginal density x-i 



Marginal density x 2 





Fig. 2 Marginal densities Ft(x) (k = 1, 2) for %i and x-i obtained from Kernel density estimation 
of random samples and from direct evaluation of Fjt(x) 



5 Numerical Examples 

5.1 Example 1 

We will use the following bivariate example with ji = (0.5,0.5)' and covari- 
ance matrix L 

1 1.2 
1.2 2 / 

as well as lower and upper truncation points a = ( — 1, —co)',b = (0.5,1)', 
i.e. X\ is doubly, while is singly truncated. The bivariate marginal den- 
sity Fq /r (x,y) is the density function itself and is shown in figure [TJ the one- 
dimensional densities F^x) (k — 1,2) in figure|2] 

The moment calculation for our example can be performed in R as 
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Monte Carlo estimator for m Monte Carlo estimator for |i 2 




2000 6000 10000 2000 6000 10000 

sample size sample size 



Fig. 3 Trace plots of the Monte Carlo estimator for fi* 



> library (tmvtnorm) 

> mu <- c(0.5, 0.5) 

> sigma <- matrix(c(l, 1.2, 1.2, 2), 2, 2) 

> a <- c(-l, -Inf) 

> b <- c(0.5, 1) 

> moments <- mtmvnorm(mean=mu, sigma=sigma, 

> lower=a, upper=b) 

and results in ji* = (—0.152, —0.388)' and covariance matrix 



/ 0.163 0.161 \ 
0.161 0.606 ) 



The trace plots in figures |3] and [4] show the evolution of a Monte Carlo esti- 
mate for the elements of the mean vector and the covariance matrix respec- 
tively for growing sample sizes. Furthermore, the 95% confidence interval 
obtained from Monte Carlo using the full sample of 10000 items is shown. All 
confidence intervals contain the true theoretical value, but Monte Carlo esti- 
mates still show substantial variation even with a sample size of 10000. Sim- 
ulation from a truncated multivariate distribution and calculating the sample 
mean or the sample covariance respectively also leads to consistent estimates 
of }i* and E*. Since the rate of convergence of the MC estimator is 0(y/n), 
one has to ensure sufficient Monte Carlo iterations in order to have a good 
approximation or to choose variance reduction techniques. 
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Monle Carlo estimator (or Cov(X(. x-i) Monte Carlo estimator for Cov{x 1: x 2 ) Monte Carlo estimator for Cov(x ? , x ? ) 




2000 S000 10000 2000 6000 10000 2000 S000 10000 

sample size sample size sample size 



Fig. 4 Trace plots of the Monte Carlo estimator for the 3 elements of E* (c n , a^ 2 = ctJj and (j| 2 ) 
5.2 Example 2 

Let// = (0,0, 0)',the covariance matrix 

/l.l 1.2 0\ 
L = 1.2 2 -0.8 
\ -0.8 3 / 



and the lower and upper truncation points a = ( — 1,-00,-00)' and = 
(0.5, 00, 00)', then the only truncated variable is X\, which is furthermore un- 
correlated with X3. 

Our formula results in ji* = c(-0.210, -0.229, 0)' and 



/ 0.174 0.190 0.0 \ 
L* = 0.190 0.898 -0.8 
\ -0.8 3.0 / 



For this special case of only k < d truncated variables {x\, . . . , X^), the remain- 
ing d — k variables (xi+i, . . . , X^) can be regressed on the truncated variables, 
and a simple formula f or the mean and covariance matrix can be given (see 
Tohnson and Kotzl Jl97l|) , p. 70). 

Let the covariance matrix L of (x\, ...,x^)he partitioned as 



(24) 
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where Vn denotes the k x k covariance matrix of {x\, . . . , x^). The mean vec- 
and the covariance matrix L* of all d variables can be computed as 

(SUiV^Vu) (25) 

and 

VVziV^Un V 22 -V 2 i(Vii 1 -Vu 1 UiiVii 1 )Vu > / ' 

where & and Un are the mean and covariance of the (x\, . . . , X/ c ) after trun- 
cation. 

The mean and standard deviation for the univariate truncated normal X\ are 

£7ll = ' 7ll + " 11 <WM-«W«i) 

Letting Un = u n and inserting £1 and Un into equations l l25t and l|26l l, one 
can verify our formula and the results for ji* and L* . However, the crux in 
using the Johnson/Kotz formula is the need to first compute the moments of 
the truncated variables [x\, . . . , x^) for k > 2. But this has been exactly the 
subject of our paper. 



6 Moment Calculation and Conditional Independence 

In this section we establish a link between our moment calcu lati on and the 
theory of graphical models ( Whittaker ( 1990), [Edwards] |l995) and Lauritzenl 
(|l996|) '). We present some properties of the inverse covariance matrix and 
show how the dependence structure of variables is affected after selection. 

Graphical modelling uses graphical representations of variables as nodes 
in a graph and dependencies among them as edges. A key concept in graphi- 
cal modelling is the conditional independence property. Two variables X and 
Y are conditional independent given a variable or a set of variables Z (nota- 
tion X _LL Y|Z), when X and Y are independent after partialling out the effect 
of Z. For conditional independent X and Y the edge between them in the 
graph is omitted and the joint density factorizes asf(x,y\z) = f(x\z)f(y\z). 

Conditional independence is equivalent to having zero elements Cl X y in 

the inverse covariance matrix CI = as well as having a zero partial co- 
variance /correlation between X and Y given the remaining variables: 

X AL Y | Rest <^=> Cl xy = <=> p xy , Rest = 

1 The formula for the truncated mean given in Johnson and Kot3 fl97ll) , p. 70 is only valid 
for a zero-mean vector or after demeaning all variables appropriately. For non-zero means fi = 
(^ 1 ,p 2 )'itwillbe(^, fe +(?i-^)V 11 1 V 12 ). 
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Both marginal independence and conditional independence between vari- 
ables simplify the computations of the truncated covariance in equation l fl6)l . 
In the presence of conditional independence of i and / given q, the terms 
~ ViqVqqVqj = vanish as they reflect the partial covariance of i and j 
given q. 

As has been shown by lMarchetti and Stan ghellini (2008]), the conditional 
independence property is preserved after selection, i.e. the inverse covariance 
matrices O and fl* before and after truncation share the same zero-elements. 
We prove that many elements of the precision matrix are invariant to trun- 
cation. For the case of k < d truncated variables, we define the set of trun- 
cated variables with T = {x\, . . . , x^}, and the remaining d — k variables as 
S = {xfc+l/ • ■ •/*£?}• We can show that the off-diagonal elements fl,-,- are in- 
variant after truncation for i 6 T U S and j G S: 

Proposition 6.1. The off-diagonal elements Cln and the diagonal elements fiy^- are 
invariant after truncation for i £ TUS and j 6 S. 

Proof. The proof is a direct application of the Johnson/Kotz formula in equa- 
tion l l2"6t in the pr evious sectio n. As a result of the formula for partitioned 
inverse matrices jGreend (|2003l) , section A.5.3), the corresponding inverse co- 
variance matrix fl of the partitioned covariance matrix £ is 



n= ( Vf/ (I + V 12 F 2 V 21 Vii 1 ) - Vii 1 V 12 F 2 
-F 2 V 21 V 1 " 1 1 F 2 



(27) 



with F 2 = (V 22 - VjxV^V^) 



Inverting the truncated covariance matrix L* in equation j26t using the 
formula for the partitioned inverse leads to the truncated precision matrix 

_ ( + V 12 F 2 V 21 V- 1 - V 12 F 2 \ 

" I -F 2 V 21 V^ F 2 J (28) 

where the fij 2 and fl 21 elements are the same as fli 2 and fl 2 i respectively. 
The same is true for the elements in Q 22 , especially the diagonal elements in 

n| 2 . □ 

Here, we prove this invariance property only for a subset of truncated 
variables. Based on our experiments we conjecture that the same is true also 
for the case of full truncation (i.e. all off-diagonal elements in O^), but we do 
not give a rigorous proof here and leave it to future research. 



6.1 Example 3 

We illustrate the invariance of the elements of the inv erse covariance m atrix 
with the famou s mathematics marks example used in IWhittakerl jl990|) and 
Id wards! dl995() . p. 49. The independence graph of the five variables ( W, V, X, 
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Y, Z) in this example takes the form of a butterfly. 

vectors (W) statistics (Z) 




mechanics (V) analysis (Y) 

Here, we have the conditional independencies (W,V) _LL (Y,Z)|X. A 
corresponding precision matrix might look like (sample data; zero-elements 
marked as "."): 



Q 



/ 1 

0.2 

0.3 
V ■ 



0.2 
1 

-0.1 



0.3 . .\ 
-0.1 . . 

1 0.4 0.5 
0.4 1 0.2 
0.5 0.2 1 / 



(29) 



After truncation in some variables (for example ( W, V, X) as —2 < W < 1, 
— 1 < V < 1, < X < 1), we apply equation lfl6l l to compute the truncated 
second moment and the inverse covariance matrix as: 



/1.88 0.2 0.3 . .\ 
0.2 3.45 -0.1 . . 
0.3 -0.1 12.67 0.4 0.5 
. 0.4 1 0.2 
V . .0.5 0.2 1 / 



(30) 



The precision matrix Q* after selection differs from Q only in the diagonal 
elements of ( W, V, X) . From Q* we can read how partial correlations between 
variables have changed due to the selection process. 

Each diagonal element Clyy of the precision mat rix is the inverse of the 
partial variance after regressing on all other variables ( Whittakerl l|l990|) .p. 143). 
Since only those diagonal elements in the precision matrix for the k < d of the 
truncated variables will change after selection, this leads to the idea to just 
compute these k elements after selection rather than the full k(k + l)/2 sym- 
metric elements in the truncated covariance matrix and applying the John- 
son/Kotz formula for the remaining d — k variables. However, the inverse 
partial variance of a scalar Y given the remaining variables X = {x\, . . . , Xj} \ 

y 



n* 



^yy ~ ^yX^XX ^Xy 



still requires the truncated covariance results derived in Section 3. 
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7 Summary 

In this paper we derived a formula for the first and second moments of 
the doubly truncated multivariate normal distribution and for their bivari- 
ate marginal density. An implementation for both formulas has been made 
available in the R statistics software as part of the tmvtnorm package. We 
linked our results to the theory of graphical models and proved an invari- 
ance property for elements of the precision matrix. Further research can deal 
with other types of truncation than we considered (e.g. elliptical). Another 
line of research can look at the moments of the doubly truncated multivari- 
ate Student-t distribution, which contains the truncated multivariate normal 
distribution as a special case. 
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